DCPT-12/47 
NSF-KITP-12-234 



Boundary conditions for scalars in Lifshitz 



Tomas Andrade'^'B and Simon F. Ross'H] 



Department of Physics, UCSB 
Santa Barbara, CA, 93106, USA 



^ Centre for Particle Theory, Department of Mathematical Sciences 
Durham University, South Road, Durham DHl 3LE, U.K. 



Abstract 



We consider the quantisation of scalar fields on a Lifshitz background, ex- 
ploring the possibility of alternative boundary conditions, allowing the slow 
falloff mode to fluctuate. We show that the scalar field with alternative bound- 
ary conditions is normalizable for a larger range of masses than in the AdS case. 
However, we then find a new instability for alternative boundary conditions, 
implying that the range of masses where alternative boundary conditions define 
a well-behaved dual theory is < < m\p + 1, analogously to the AdS 

case. The instability is of a novel type, with modes of arbitrarily large momen- 
tum which grow exponentially in time; it is therefore essentially a UV effect, 
and implies that the dual field theory is simply not defined where it appears. 
We discuss the interpretation in the dual field theory, and give a proposed lower 
bound on the dimension of scalar operators. 

1 Introduction 

The holographic description of field theories with anisotropic scaling symmetry presents 
an interesting extension of AdS/CFT, which may have valuable applications in con- 
densed matter theory. The simplest example of such a dual description is the Lifshitz 
metric originally constructed in pLj. The geometry is 



where represents the overall curvature scale, and the spacetime has d + 1 di- 
mensions, so there are ds = d — 1 spatial dimensions x. The asymptotically Lif- 
shitz solutions of a bulk gravity theory are conjectured to provide a dual holographic 
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description of a non-relativistic field theory, with the anisotropic scaling symmetry 
t — X^t, X* — Xx\ This duality is interesting both for its potential application 
in condensed matter, and as an extension of our understanding of holography and 
the relation between field theory and spacetime descriptions. The holographic dic- 
tionary, relating bulk spacetime quantities to field theory observables, is now fairly 
well-developed PQI21E]. As in AdS/CFT, this identifies the leading asymptotic falloff 
of bulk fields with sources in the dual field theory. 

An interesting early observation in AdS holography p[j was that there can be 
more than one conformal field theory associated to a given bulk theory, depending 
on the choice of boundary conditions. As first noted in [5] for scalar fields, for some 
parameter values it is possible to introduce an alternative quantisation in the bulk 
spacetime, where a mode which is subleading in the asymptotic expansion of the field 
is fixed and the leading piece is allowed to vary. This alternative quantisation leads 
to a second conformal field theory dual to the same spacetime, with different operator 
dimensions for the operators dual to the bulk fields whose boundary conditions have 
been changed. One can also consider mixed boundary conditions, which are dual to 
renormalisation group flows interpolating between the two conformal field theories, 
generated by a double-trace deformation of the field theory. For a scalar field, the 
alternative quantisation is possible when the mass of the field is in the range m%p < 
m? < m\p + 1, where m\p = — ^ is the Breitenlohner-Freedman bound [5] for d 
boundary dimensions. 

The importance of this possibility for holography was realised in part because 
the dimension of the operator dual to the scalar with standard boundary conditions, 

A_|_ = I + \J ^ + '^^j was always strictly greater than the unitarity bound, A > 

I — 1. With the alternative boundary condition, the scalar is dual to an operator 

with dimension A_ = ^ — ^J~^^^\-rn^, which precisely saturates the unitarity bound 

when = m%p + 1. 

This analysis of alternative boundary conditions in AdS deepened our understand- 
ing of the correspondence, through an improved bulk understanding of unitarity and 
an understanding of the relation of double-trace and more general deformations of 
the field theory and boundary conditions in the bulk [HI E] . The early work on scalar 
fields was subsequently extended to consider vector fields in [51 IH], and for metric 
perturbations in [TO] . 

In the bulk spacetime, the restriction to the range < rn^p + 1 comes from 
the fact that the slow fall-off mode is only normalizable with respect to the usual 
Klein-Gordon norm for masses in this range. Similar restrictions arise for a vector 
field [mil]. A deeper understanding of this restriction was obtained in [TU], where it 
was observed that one can define an alternative norm by adding boundary terms to 
the Klein-Gordon current, such that the solution with Neumann boundary conditions 
is always normalizable. However, this norm is not positive for > -|- 1 for 
generic boundary metrics [10], and for the flat boundary metric, where the norm is 
positive, an IR divergence appears enforcing the unitarity bound fH]. 

In this paper, we aim to address the same issues in Lifshitz. We will focus primarily 
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on investigating when alternative boundary conditions are possible from the bulk 
spacetime perspective. In the present paper we do this calculation for a scalar field 
in a fixed Lifshitz background; a companion paper to follow will consider linearised 
fluctuations in a theory including dynamical gravity. Our intention is to use this 
study to shed further light on the duality. In particular, the calculation with a scalar 
field will lead to a new prediction for a bound on operator dimensions in Lifshitz field 
theories (or at least those with gravitational duals). 

After some preliminary discussion of the equation of motion in section [2| we first 
investigate the normalizability for a scalar field on a Lifshitz background in section 
[3| We find that the slow fall-off mode is normalizable, so that alternative boundary 
conditions are possible, for a larger range of masses than in the AdS case: 

m%p < < m%p + , (2) 

where m\p = —\{z + dsY and ds is the number of spatial dimensions. The lower 
bound m\p is the analogue of the Breitenlohner-Freedman bound in this case. The 
immediate source of this wider range is that the measure in the Klein- Gordon inner 
product includes a factor of r~^, so slower falloffs become normalizable as z increases. 
More deeply, this can be related to the perspective of [10] by observing that as z 
increases the dimension of boundary counterterms involving time derivatives of the 
boundary data increase, and the inner product needs to be modified by adding bound- 
ary terms only when the operator dimension is small enough for these terms to be 
relevant. In the field theory, the wider mass range ([2]) corresponds to a weaker restric- 
tion on the possible dimensions of scalar operators. Thus, the spacetime calculation 
predicts a reduced lower bound for the dimension of scalar operators in the field 
theory. 

However, in section |4], we show that the scalar field on a Lifshitz background with 
alternative boundary conditions can be unstable even when m? > m?Bp. In the AdS 
case, an instability for conformally-invariant boundary conditions is kinematically 
forbidden: conformal symmetry implies that the dispersion relation is a; = ±fc. How- 
ever for Lifshitz boundary conditions, the scaling symmetry only fixes the dispersion 
relation to be a; = ak^ for some dimensionless parameter a. If there is a family of 
modes where a has a positive imaginary part, they represent an instability. 

We will first consider the question analytically for 2 = 2, where the solution of 
the scalar equation can be written in terms of confluent hypergeometric functions. 
We show that the theory with Dirichlet boundary conditions is stable, but that the 
theory with Neumann boundary conditions is unstable if 

> m|p + 1. (3) 

We note that the instability is a UV effect, as if there is any instability then there will 
be exponentially growing modes for all values of /c, and the dominant instability is 
associated with arbitrarily large values of k. Thus, we would expect this instability to 
appear for any asymptotically Lifshitz spacetime when we take Neumann boundary 
conditions for the scalar. We confirm this expectation analytically by considering the 
scalar field on the z = 2 Lifshitz black hole solution introduced in [13], where the 
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wave equation can be solved in terms of hypergeometric functions. This makes it 
clear that the instability we are seeing here is not associated with the singularity in 
the interior of the spacetime for the pure Lifshitz solution [U [HI [15] ; this is a new UV 
pathology in these solutions. We complete our analytic discussion by giving a general 
argument for stability of black hole solutions with Dirichlet boundary conditions, 
following arguments previously given in [T3] . 

We then extend our discussion to arbitrary values of z by solving the wave equation 
numerically. We do numerical analysis both for the pure Lifshitz spacetime and for a 
spacetime with an IR cutoff, which is easier to treat numerically, using both a spectral 
method and shooting. We verify that the IR cutoff does not alter the instability for 
high momentum. We find that the theory with Neumann boundary conditions is 
unstable if > + 1 for all z. 

In the field theory dual to the alternative quantisation, the dimension is given by 
A_ = \(ds ^ z) — \ \/ [ds^ zY + 4m^, so w? < m^^p + 1 corresponds to 

A>^(4 + ^)-l. (4) 

This generalises the result A > | — 1 in the relativistic case (where d = + 1 the 
number of spacetime dimensions) . In the Lifshitz field theory, there is no independent 
derivation of such a bound; in the absence of the usual state-operator map for such 
non-relativistic theories, there is no direct analogue of the usual argument for the 
unitarity bound. However, it is certainly surprising that the bound is so high; if 
we consider a free scalar field theory with kinetic term {dt4>)^, Lifshitz scaling would 
require the scalar to have dimension A = ^{ds + z) — z. This matches the bound that 
would be obtained from normalizability considerations alone, but the instability we 
find suggests that for interacting Lifshitz field theories the bound on the dimension is 
higher. We will discuss the field theory interpretation a little more in the conclusions 
in section [5l 

Note added: While this paper was in preparation [T7] appeared, which has con- 
siderable overlap with our analysis of normalizability in section [3j 



2 Scalar wave equation 

We consider a scalar field satisfying the Klein-Gordon equation 00 — = in the 
Lifshitz spacetime ([T]) with < 0, and we neglect back-reaction on the spacetime 
metricQ The wave equation in the Lifshitz geometry ([T]) is 

r^-^-^^9,(r"+^^+^9,0) - {r-^'d^ - r'^d^ + m^)^ = 0. (5) 

As in AdS, the derivatives along the boundary direction have a subleading effect at 
large r, and the asymptotic behaviour of the solutions is 

0~ 0+r-^+ +0_r-^-, (6) 
^Extending the analysis to include back-reaction might be an interesting project for the future. 
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where 

Ai = ^{z + ds)±^m^ + ^{z + dsy. (7) 

The analogue of the Breitenlohner-Freedman bound [5] for Lifshitz spacetimes is thus 
rri^p = — \{z + dsY. We will also write = m'^p + , so A-t = \{z \ d^±i z/, with 
z/ > by convention. The usual boundary condition is to fix the slow fall-off mode 0_ 
and let 0+ fluctuate. We will generally refer to this as a Dirichlet boundary condition, 
and to the converse condition of fixing 0+ and letting 0_ fluctuate as a Neumann 
boundary condition, although this terminology is really only valid for m = where 
A_ =0. 

To solve the equation (|5]) explicitly, we use a boundary plane wave basis, writing 

= e-*'^*+^^->(r). (8) 

The wave equation then becomes 

r^-^-'^=9,(r"+'^=+^9^V^) - {-r-^'u? + r-'^k^ + m^)^ = 0. (9) 

The dg dependence here can be simplified by writing '^(r) = r~^^x(r); then 

rdrirdrX) - {-r'^'co^ + r'^k"^ + i^^)x = 0- (10) 

The asymptotic behaviour of x is then easily seen to be x ~ ^^'^ as r — )■ oo, with the 
plus (minus) sign corresponding to Neumann (Dirichlet) boundary conditions, and 
~ e^*^ as r — )■ 0. For real u both behaviours at r — )■ are regular. For complex 
u one grows exponentially and the other decays; we select the exponentially damped 
mode. 

We can also note that either u or k can be absorbed by a redefinition of r, so the 
equation can be rewritten as 

rdrirdrX) - + - r"'^^^ X = 0. (11) 

This makes manifest the fact that the physics on the pure Lifshitz background can 
depend only on the dimensionless combination u/k^. 



For generic u and z, (10) has no solution in terms of known special functions. For 
w = 0, we can solve it in terms of Bessel functions, but for both the normalizability 
and instability discussions our interest is in solutions with non-zero u. More helpfully, 
for z = 2, the equation reduces to a confluent hypergeometric equation; this was 
previously analysed in [1], and we will use this solution in studying the instability in 
section HI 



3 Normalizability for scalars 

In this section, we consider the normalizability of the probe scalar field. We first give a 
simple consideration of normalizability with respect to the normal Klein-Gordon inner 
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product, showing that the range of masses for which alternative boundary conditions 
are possible is enlarged as we increase z. We then do a more detailed analysis of 
the inner product and counterterms, following [101 12] closely. This allows us to 
understand the result better from the spacetime point of view, seeing the relation 
to kinetic counterterms in the action. We verify that inside our mass range, the 
standard Klein-Gordon inner product without any explicit boundary contributions is 
an appropriate inner product; in particular it is finite and positive definite for real uj. 
We assume we use the standard Klein-Gordon inner product. 



d'''xdr^n^{<^\d^<^2 - 02^^01), (12) 



where S is a spacelike surface, which we will take to be a surface of constant t. The 
wave equation (|9| can be written as a Sturm-Liouville (SL) problem with eigenvalue 
\ = ikP' for the operator 



-|-(Mr)|-)+g(r) 



(13) 

with p = r^'^s^^^'^) ^ uj = ^ds-z-i _ ^.{ds+z-i) _|_ ^-2^2^ _ rjj-^g inner product 



(12) then becomes 



(01,02) = {2nr^S^'^\h - fc,)e-K-^)* ^^^+^^^ (^i,^2)5L, (14) 
where {■,-)sl is the corresponding SL inner product 



oo 



(^i,^2)5L= / r'^=-^-Vt^2C^r-. (15) 
Jo 

With the standard Dirichlet boundary conditions 0_ =0, the fields fall off as 
~ r~^+ , and the large r behaviour of the integral is J°° r^'^'^~^^~^dr, so the fast fall 
off modes are normalizable for any u. If we consider instead the Neumann bound- 
ary condition 0+ = 0, then ~ r~^- , and the large r behaviour of the integral 
is J°° r'^'^~'^^~^dr, so the slow fall off modes are normalizable with respect to this 
standard inner product if u < z, that is if 

rn^BF < < "^BF + -2^- (16) 

Increasing z thus increases the mass range for which the Neumann boundary condi- 
tions are allowed. In this range, we could also consider mixed boundary conditions; 
the flux through infinity vanishes, so that the inner product is conserved, for any 
hnear boundary condition 0+ = /0_ for real /. 

The theory with Neumann boundary conditions is dual to a field theory with 
Lifshitz scaling with an operator of dimension 

A_ = ^{z + ds)-u. (17) 
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The more general mixed boundary conditions are dual to non-scale invariant theories, 
which should interpolate between the Neumann boundary conditions in the UV and 
Dirichlet boundary conditions in the IR. We note that for z > dg, A_ could formally 
take negative values. We will find an instability in the next section before we reach 
such values, however, so we will confine ourselves to considering u < \{z + dg), where 
a linearised analysis of the asymptotics is possible. 



3.1 Inner product and counterterms 

To gain a better understanding of the origins of the extended region of normalizability, 
we should consider a more careful analysis of the inner product, following [TnilI2]- The 
argument above assumed that we could use the standard Klein-Gordon inner product. 
However, pTO] made it clear that the presence of kinetic boundary counterterms in the 
action for fields in AdS implies that we generally need to add corresponding boundary 
terms to the inner product to ensure that it is conserved. The correct inner product 
is always finite, as the boundary terms cancel any divergence from the bulk, but these 
terms may spoil positivity. The correct question to ask is then when we need to add 
boundary contributions to the Klein-Gordon inner product, and to check that the 
inner product is positive and conserved. We will see below that < m\p + z'^ 
is precisely where no explicit boundary contribution is required to have a conserved 
inner product. 

The terms (j)± in ^ can be expanded in a double power series in r~^, r~^^. 



where the terms which involve powers of are local functions of cf)'^^'' and its spatial 



derivatives, while the terms which involve powers of r are functions of temporal 
derivatives of For example, 

while 

\i u < z^ the terms involving temporal derivatives of (j)^^'^ are subleading compared to 
0*^^") in the asymptotic expansion of 0: 

^ ^-|(-+rf.)+-(0(o) + ^-2^(1) ^ ^-4^(2) ^ ^ _ ^ ^ r-2>(-) + r-2^0(^) + ...). (21) 

To obtain a finite on-shell action, one needs to add counterterms to the bare action 
to cancel divergences associated to the terms 0*^°-' and 0*-*^ for i < u. For example. 



^We don't write the subleading terms involving derivatives of 4>''^^ explicitly because they don't 
enter the calculation. 
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if 1 < z/ < 2, there will be divergences in the on-shell action involving both 0*^°^ and 
0(1), 



' bare 



(f^xdt^/^n^'(t)d^(t) (22) 



r=ro 



- I rf'^=xrftr2'^[A_0(o)2 + (A_ + 2)0(°V^'Vo-2 + A_ 

2 J r=ro 



The divergences are cancelled by adding appropriate counterterms. The leading di- 
vergence is cancelled by a 0^ counterterm, and the (p^^^cp^^^ divergence can be cancelled 
by a {di(j)Y counterterm. 

In [lOj, it was shown that derivative counterterms in the action naturally lead 
to boundary contributions to the symplectic structure. At the level of the inner 
product, we can understand these terms as being required to ensure conservation 
of the inner product: we add boundary terms to the Klein-Gordon current so that 
the flux through the boundary at large r vanishes. The key difference in our case is 
that the counterterm in the action only involves spatial derivatives. As a result, the 
usual Klein-Gordon current will be conserved despite the appearance of a (j)'^^^ term 
in the flux through the boundary. The flux of the Klein-Gordon current through the 
boundary is 



F = I [ d''^xdt^n^{<j)ld^<f)2 - 02^^01) 

J r=ro 



(23) 



J r=ro 

If we work in the plane-wave basis, the divergent term is 

d'^^xdtrl''-^ ^1 ~ e^('-^-'-^)*--(^'^-^^^)-">i°)'^(r)4°Hr), (24) 

and if we integrate over the region between two surfaces t = constant, the integral 
over the spatial directions will introduce an overall momentum delta-function, so that 
this divergent term vanishes]^ As claimed earlier, the finite piece will vanish for any 
mixed boundary condition 0+ = /0_ with a real coefficient /. Thus, for u < 2, the 
Klein-Gordon inner product is conserved for any such mixed boundary conditions; 
in particular it is conserved for both Dirichlet and Neumann boundary conditions, 
without adding an explicit boundary term to the inner product. 

This argument can easily be extended to general u < z. The key point is that 
the subleading terms appearing in the asymptotic expansion for will all involve 
only spatial derivatives of 0(''), so they make vanishing contributions to the total flux 
through the boundary at infinity between two surfaces t = constant. 

Once we know that the inner product is conserved, it is easy to see that it is 
orthogonal in the plane wave basis = e~*'^*"'"*'^''^'0(r). We saw already in (14) that 



■^One can add a local counterterm to cancel this divergence; this counterterm will not affect the 
value of the inner product. 
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the spatial integral makes the inner product vanish if ki ^ Now the fact that it 
is independent of the spatial slice t = constant we choose to evaluate it on implies it 
must vanish if cui 7^ U2 for real u, as otherwise it would be time dependent. We are 
therefore left with the inner product of the plane wave modes with themselves, 

POO 

{(pi, (pi) = {27ry^V ol{x)uJi drr'^^-'-'^\ipi\'^, (25) 

Jo 

where Vol{x) is the spatial volume. This is manifestly positive if > 0. 

We can therefore understand the condition u < z as arising from requiring that 
the free data 0'-^^ appears in the asymptotic expansion before the first term which 
involves time derivatives of 0*^°^. That is, it is precisely because we do not require 
counterterms in the action involving time derivatives that the standard Klein-Gordon 
norm remains appropriate for more general boundary conditions. 



4 Instability for Neumann boundary conditions 

We would now like to consider the spectrum for the different possible boundary con- 
ditions, to check if there are any instabilities, looking for regular solutions which grow 
exponentially in time. From the point of view of the dual field theory, instabilities 
appear as poles in the two-point function in the upper half frequency planej^ 

For the AdS case, the symmetries imply that the two-point function is a function 
of the Lorentz invariant u"^ — k'^ . Instabilities can then occur only in the case of 
mixed boundary conditions, where they correspond to tachyonic poles in the two- 
point function with uj"^ — k"^ = —fnl^y^ For the conformally invariant pure Dirichlet 
or Neumann boundary conditions, by contrast, no such instability is possible, as there 
is no scale to provide a value for m^^^. 

In the Lifshitz case however, the symmetry is less restrictive, and the two-point 
function for general boundary conditions can depend separately on u and k. For con- 
formally invariant boundary conditions, the two-point function (up to overall scaling) 
must be a function of the invariant uj/k^, but this still admits the possibility of in- 
stability, if the two-point function has a pole at 

^ = a, Im a > 0. (26) 

k^ 

We will see below that such unstable modes appear for Neumann boundary conditions 
when V > 1. The existence of such instabilities in the conformally invariant case is 
perhaps surprising. Moreover, as a result of the scale invariance, the instability has no 
associated timescale. Unlike the relativistic case, when an instability occurs, there are 
exponentially growing modes for all momenta, and the modes of high momenta have 

^ Since our situation is time-translation invariant, poles at complex will appear in complex 
conjugate pairs. 

^The mixed boundary conditions correspond to the field theory deformed by a multi-trace opera- 
tor, and a bulk instability can be interpreted as this deformation making the field theory Hamiltonian 
unbounded from below [HI HH] . 
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arbitrarily high growth rates. Thus, this instabihty indicates the complete breakdown 
of the expansion around the putative background, which is not valid even in an open 
neighbourhood in time. We therefore interpret this instability as saying that no dual 
field theory exists for the boundary conditions where the instability is present. 

Because there are unstable modes of arbitrarily high momentum, this is essentially 
a UV effect in the field theory. That is, this instability will affect not just the pure 
Lifshitz spacetime, but any solution which asymptotically approaches this solution 
with the Neumann boundary conditions. 

In particular, this implies that the theories with mixed boundary conditions will 
also be unstable, as these approach the Neumann boundary condition in the UV. That 
is, the field theories dual to the mixed boundary conditions are relevant deformations 
of the scale-invariant theory dual to the Neumann boundary conditions, so the non- 
existence of this UV theory implies that the theories dual to the mixed boundary 
conditions will also not be well-defined. 



4.1 Analytic calculation for z = 2 

We first consider the case z = 2, where it is possible to write the solution explicitly 
in terms of confluent hypergeometric functions, as discussed in pj. For z = 2, the 
generic solution of (|9|) is 



^{u) = e'''''^^^[^+u-^+\F^{a, b, -iu^u) + ^^u-^-\Fi{a -b+l,2-b, -iu^oj)], 

(27) 

where u = 1/r, 

1 k"^ 

a= -(l + z/)+i — , b=l + v, (28) 
2^ ' Auj 

and iFi(q;,7, z) is the confluent hypergeometric function, whose series expansion is 
iFiia, 7, 2;) = \ + ^z + . . .. The first term in (27) corresponds to the fast falloff mode 
at infinity, and the second term to the slow falloff. 

For complex the solution which is regular at r = is the one which is ex- 
ponentially damped. The solution which is regular for Im a; > can be written 
as 

^ ^ g.„V2^^+.^ _,^„2^ (29) 



where u = 1/r, the constants a, b are given in (28 ), and U (a, 6, z) is Tricomi's confluent 
hypergeometric function, which is given in terms of the confluent hypergeometric 
function by 

^' ^) = f^r^^TT^i^i(«' ^' ^) + ^pr/^ ^'AFi(a - 6 + 1, 2 - 6, z). (30) 
1 (a — + Ij 1 (aj 

We then want to ask for values of u such that this also satisfies the boundary condition 
at M = 0. Using (30), the solution near -u = is 



^ = -"^-^^#^(1 + •••)+ u-^-^ n^^'K.M + •••), (31) 
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and it will satisfy Dirichlet (Neumann) boundary conditions if the first (second) term 
vanishes, which is when the corresponding Gamma functions have a pole. This implies 
that the frequencies must satisfy 

" = -' 2(l + 2l±.) = 0.1.2.3.- (32) 

where the plus (minus) sign corresponds to Dirichlet (Neumann) boundary conditions. 
With Dirichlet boundary conditions this condition cannot be satisfied with Im w > 
for any n. There are thus no solutions which satisfy Dirichlet boundary conditions 
and are regular on the horizon for Im 7^ 0. But for Neumann boundary conditions 



and z/ > 1, there are solutions of (32) with Im a; > 0. These are therefore unstable 
modes which satisfy the boundary conditions both at infinity and the horizon. We 
see that the number of unstable modes increases as v crosses integer values. 

Note that u is pure imaginary; unlike the usual quasi-normal modes in the rela- 
tivistic case which have decaying oscillations, the unstable mode simply grows expo- 
nentially with no oscillatory component. We can thus also view this as the analytic 
continuation of a mode of real frequency on the Euclidean Lifshitz solution. 

4.1.1 Black hole for 2 = 2 

As we argued above, this instability should appear in any asymptotically Lifshitz 
spacetime. One simple case where we can check this explicitly is for the black hole 
considered in |13j^ The geometry is 

ds^ = -r^ (l - dt^ + — + r'^dx^. (33) 



This metric can be obtained as a solution of a higher-derivative theory of gravity 
[13] . but for our present purposes the point is that it provides an example of an 
asymptotically Lifshitz spacetime where the wave equation can be solved analytically. 
For this metric the wave equation in the boundary plane wave basis becomes 

r^'-'dr{r'^^^\r' - r^)9.V^) - (-^ + ^ + ^ = ^- ^ 

Introducing a radial variable y = ^2 ^ , this becomes a hypergeometric equation |20] . 
The solution which is well-behaved at the black hole horizon y = for Im w > is 

^ = t/"^(l -y)j('^+')+5-F(a,/3,7;2/), (35) 
where F[a,P,'j] y) is the hypergeometric function and 



^Thc ds = 2 version of this metric was previously found in 20J as a solution of Einstein gravity 
coupled to various matter fields. The scalar quasi-normal mode spectrum was discussed also in |20] 
for particular values of the couplings, and then generalized in |21) . 
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Analytically continuing the hypergeometric function to y = 1 gives 



r(7)r(7-«-/3) 
r(7-«)r(7-/3) 



F{a,(3,a + (3- J + 1,1- y) 



(37) 



+(1 _ ^)7~»-/3 r(7)r(a + ^ 7) ^^^ _ ^ _ ^ _ « _ ;3 + 1, 1 _ 

r(a)r(/:/) 

Thus, the solution will satisfy Dirichlet boundary conditions at infinity if the second 
term vanishes, which is to say if there are poles in the Gamma functions in the 
denominator, a = —n or /3 = —n for some integer n. Rearranging and squaring gives 



2iu{l + iy + 2n) = k'^ + rl(^{l + iy + 2nf - ^ . 



(3J 



We see that for large k, the term proportional to ri becomes negligible, and the 



spectrum reduces to the one found in (32), as expected. There is thus no instability 



at large k; however it is interesting to note that for ds > 2, the black hole can produce 
instabilities for small if z/ is small enough. The black hole will have at least one 
exponentially growing mode if (1 + z/)^ — ^ < O/] 



The solution (35) will satisfy Neumann boundary conditions if the first term in 
7) vanishes, that is if 7 — a = —n or 7 — /3 = —n for some integer n. As in the 
pure Lifshitz spacetime this simply corresponds to — )■ —v in the Dirichlet analysis 
above. The spectrum of such modes is given by 



2iuj{l -u + 2n) = P + rl[{l-u + 2nf - 



dl 



(39) 



Again, this reduces to (32) at large k. We see explicitly that the unstable modes at 



large k are present for the black hole solution as well. 



4.2 Stability for Dirichlet boundary conditions for general z 

For general z, we cannot solve the wave equation analytically. However, for Dirichlet 
boundary conditions, we can give a general argument that no such instability can 
occur. This argument does not consider the pure Lifshitz spacetime directly; instead 
it is convenient to regulate the IR behaviour by considering an asymptotically Lifshitz 
black hole solution and using an argument for stability originally developed for AdS 
black holes in [16] and extended to the Lifshitz case in [13]. The absence of an 
instability for large k on the black hole spacetime will then imply the absence of any 
instability on the pure Lifshitz spacetime. 

We consider the black hole in ingoing Eddington-Finkelstein coordinates, where 
we assume the metric is of the form 

ds^ = -r^^f(r)dv^ + 2r'-^dvdr + r^dx"^, (40) 

^This instability was missed in I13j because they only considered Dirichlet boundary conditions 
with > 0, corresponding to v > ^{ds + 2). 
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where /(r) is some emblackening function with f{r^) = 0, f'{r) > for r > r_|_, and 
/(r) — 7- 1 as r — 7- oo. As we are introducing the black hole simply as a convenient IR 
regulator, we do not specify a precise form for /(r) or consider the equations of motion 
for the metric. We will take the scalar field to be of the form = e~*'^'"+*'^'^r~^i?(r). 
Then the wave equation becomes 

{r'+^f{r)Rj - 2iujB! - V{r)R = 0, (41) 

where 

V{r) = er--' ^^^^^^(1 - /(r)) + |r/] r^-\ (42) 



To argue for stability, we want to see that Im a; < 0. By multiplying (41 ) by R* 
and integrating, we can see that 



POO 

/ dr[R*{r'+^f{r)Ry - V{r)\R\^ - 2iuR*R'] = 0. 

J r-i- 



(43) 



Thus, integrating by parts gives that 



X=l dr[r'+^f{r)\R'\^ + V{r)\R\^] = -2tuj drR*R' + r'+^f{r)R*R'\^ . (AA) 



Now f{r+) = ensures that the boundary term at the horizon vanishes (this is why 



it is useful to introduce a black hole in the bulk), and at r — oo, the solution of (41) 
asymptotically behaves as i? ~ r~5='='^. The lower sign is the fast fall-off mode which 
is retained for Dirichlet boundary conditions, so for Dirichlet boundary conditions, 
the boundary term at infinity also vanishes. Thus, 

POO 

X = -2iu / drR*R'. (45) 



Since X is real, taking the imaginary part of this equation gives (cj — cu) drR*R' 
u\R\l_^. , so 



Im a; + 
So if X > 0, Im u < and the mode is stable. 

For sufficiently large u, V{r) is positive everywhere, and X is obviously positive. 
However, we wish to give an argument for all positive u. This is possible because we 
are concerned only with instabilities which continue to appear for arbitrarily large k: 
only such instabilities can correspond to instabilities of the pure Lifshitz spacetime. 
We can therefore focus on the region at large r in the integral defining X; if there is 
a problem at large k it must come from this part where the k"^ contribution in V is 
suppressed relative to the other terms. The large r part of X is approximately 



/oo 



.2 

4 



dr[iu + + iu' - -)]r- 

) 

dr[2u'^ + uz]r-^''-\ (47) 



v2 
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where we used the falloff R ~ r~2~'^ for Dirichlet boundary conditions again in the 
second hne. Thus this large r contribution is positive and bounded, and at least 
for sufficiently large k, X > and there will be no instability. We do not exclude 
the possibility of an instability at finite k -we saw in the previous section that such 
instabilities can in fact occur - but the absence of instability at arbitrarily large 
k implies the pure Lifshitz solution with Dirichlet boundary conditions cannot be 
unstable. 



4.3 Numerical analysis for Neumann boundary conditions 
for general z 

Having seen that the theory with Dirichlet boundary conditions is stable, we now 
turn to a numerical analysis to determine when the theory with Neumann boundary 
conditions is unstable. The numerical determination of the spectrum is conceptually 
the same as the analytic calculation for z = 2: we choose a solution which is regular in 
the interior and look for values of u which satisfy the asymptotic Neumann boundary 
conditions. This defines an eigenvalue problem for u, which we can solve either by 
spectral methods or by shooting. 

For integer z, we can apply spectral methods to solve the problem in the pure 
Lifshitz geometry, as the solution is analytic even near r = after appropriate redefi- 
nitions. For non-integer values of z, we need to modify the problem, as the numerical 
analysis for pure Lifshitz is difficult because of the behaviour near the horizon. For- 
tunately, we have seen that the instability that we are interested in is essentially a 
UV effect in the field theory, so we can modify the geometry in the interior of the 
spacetime without affecting the instabihty. We could consider a black hole solution, 
as we did previously, but for numerical studies it is more convenient to simply cut 
off the spacetime at a hard wall. We introduce a radial cutoff, rescaling the radial 
coordinate so that the cutoff is at r = 1 and consider the region r > 1. For simplicity, 
we focus on Dirichlet or Neumann boundary conditions at the wall; these two choices 
lead to similar results, and give the same region of instability in u. 

We have done the numerical analysis for the theory with an IR wall using both 
spectral methods and shooting. As the wall breaks the radial scaling symmetry, we 
have to consider varying u and k separately. We are interested in solutions at large 
k, giving instabilities of the pure Lifshitz spacetime. We solve the eigenvalue problem 
for u"^ for fixed k, and check that k^/u is invariant as we vary k, for relatively large 
k. Failure to satisfy this property is interpreted as a sign that the numerics are 
inaccurate. 

As an illustration, the values obtained by shooting for z = 2 are plotted in figure 
[l|p| We see that k'^/u does indeed become constant at large k, and that for large values 



of k, the relation (32 ) is accurately reproduced. In particular, as we decrease u towards 
1/ = 1, the approximately invariant quantity k'^/u goes to zero. Representative plots 

®Note that even though the radial equation can be solved analytically in this case, the spectrum 
for the hard wall must be determined numerically since the boundary condition at the wall involves 
(roughly speaking) determining zeros of the confluent hypergeometric function. 
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Figure 1: Unstable modes for z = 2 for Neumann boundary conditions at the IR wall, 
as a function of k for various values of u. The lines of data correspond to u = 1.1, 
u = 1.15, u = 1.20 and u = 1.25, with u increasing as we move up the plot. We 



observe that for sufficiently large k, the relation (32) is reproduced. 



for different values of z using spectral methods are given in figure [2j Plots obtained 
by shooting are given in figure [3| We see clearly that there is an instability for i/ > 1. 

A significant difference between the numerical results and the behaviour of the 
analytic solution for z = 2 is the behaviour when u > z for z < 2. We find that as 
u ^ z, the invariant u/k^ calculated for large k goes to zero. This can clearly be seen 
in the plots generated by shooting in figure |3j k^ /co diverges as u ^ z (similar results 
are obtained by spectral methods for z = 1.5). The general behaviour we find is that 
there are no complex frequency modes at large k for u > z, the only mode solutions 
are those with real frequencies. There are however complex frequency unstable modes 
at small k in the IR wall geometry. 

In the limit as z — )■ 1, all the frequencies smoothly approach those in the IR wall 
cutoff version of AdS; in particular, there is a tachyonic instability in the IR cutoff 
AdS with Neumann boundary conditions whose scale is set by the location of the wall, 
and this is reproduced by the behaviour of the complex frequency unstable solutions 
at small k for u > z. Some representative large k values are plotted in figure 4. 

This explains the relation to the relativistic case: if we were to consider varying z 
at fixed u, the UV instability will disappear at u = z, and in the region near z = 1 the 
qualitative behaviour is similar to the AdS case. Note that while the UV instability 
disappears, since u > z, we are violating the normalizability bound and there will be 
pathologies for Neumann boundary conditions along the lines of the analysis in |12] . 
In particular, there will be ghosts in the spectrum in the IR wall geometry. 



5 Discussion 

In this paper, we have considered the possibility of alternative boundary conditions 
for scalar fields in the Lifshitz spacetime. We found that the usual normalizability 
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Figure 2: Scale invariant combination ik^ /uj for the first unstable mode as a function 
of V. The top (black) data is 2; = 3/2, the dotted curve (green) is the analytic result 
for z = 2, the lower data are z = 5/2 (red) and the lowest line of data (blue) is z = 3. 
The numerics become inaccurate for large problems appearing roughly for z/ > 1.4 
for z > 2. These are identified by varying the size of the grid and noticing that the 
eigenvalues change considerably. 



condition implied that Neumann boundary conditions are possible for a larger range 
of masses than in the AdS case. However, we found that there is a new instability 
for scalar fields with Neumann boundary conditions on Lifshitz, which implies that 
there is a well-defined dual field theory only for m? < m\p + 1. This instability is 
a UV effect; we saw explicitly that it is not tied to the presence of a singularity at 
r — )■ 00 in the Lifshitz geometry as the same instability appears in an asymptotically 
Lifshitz black hole and in a spacetime with an IR hard wall cutoff. 



This led to the proposed lower bound (17) on the dimension of scalar operators 
in strongly-coupled field theories with Lifshitz scaling symmetry. It would be very 
interesting to understand this bound better from the dual field theory perspective. 
The bound is non-trivial: the field theory of a free scalar field, with the kinetic term 
J dtcfi" x{dt4>Y 1 '^ill be invariant under a Lifshitz scaling with dynamical exponent z 
if the scalar field has dimension = }^{ds — z). This free scalar value violates our 
bound (p!7|)p| so if the proposed bound is to be valid, the interacting nature of the 
theory must play a central role. 

The bound for z = 1 corresponds to the unitarity bound in the field theory. 
Could there be a similarly general derivation for z > 11 The usual argument for 
the unitarity bound is based on considering the radial Hamiltonian and relating the 
dimensions of operators to the norms of states. With Lifshitz symmetry there are 
no special conformal transformations, so we cannot construct such a state-operator 



^Although intriguingly it precisely saturates the bound at which normaUzabihty fails to apply, 
when IV? — m^p + in the bulk. 
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Figure 3: We plot ik^ /uj for large as a function of v ioi z = 1.3 (top), z = 1.4 
(center) and z = 1.5 (bottom), using shooting with an IR wall. Note that ik^ /uj 
diverges as p approaches z from below. 



mapping. However, a different derivation of the unitarity bound which uses the optical 
theorem, and does not rely directly on a state-operator mapping, was given in [22] . 
This argument could in principle be extended to the Lifshitz case to obtain a bound 
directly from the field theory point of view, as was recently done for time-dependent 
cases in [23]. Unfortunately, to do this in practice seems technically challenging. The 
argument was based on relating the imaginary part of the two-point function of an 
operator O to the optical theorem for some external field coupled to O. It is possible 
to calculate the two-point function in momentum space for the case z = 2. But to 
apply the argument of [22], we need to control the overall normalisation of this two- 
point function, as function of m? and z. In [22], the overall normalisation is fixed by 
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Figure 4: We plot the frequencies of the first two modes for v = 1.3, k = 3.0 as a 
function of z with z < u. All the frequencies are real and they approach the relativistic 
values as 2; — )■ 1. 
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requiring that {O^ {x)0{y)) is positive definite at spacelike separation. Unfortunately, 
it is not possible to analytically compute the integral required to evaluate the two- 
point function obtained for z = 2 in position space. Further exploration of the field 
theory is left as an open problem for the future. 

It would also be interesting to explore these questions of alternative boundary 
conditions and instability for the full theory with dynamical gravity. We intend to 
explore these issues for linearised fluctuations of the metric and matter fields in the 
context of the massive vector model of [2l] in the future. 
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